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In this short note we discuss the relation between the so-called Off-Diagonal-Long-Range-Order 
in many-body interacting quantum systems introduced by C. N. Yang in Rev. Mod. Phys. 34, 694 
(1962) and entanglement. We argue that there is a direct relation between these two concepts. 



We all know that entanglement is useful in quantum 
information processing. In principle, it is possible to built 
a machine (quantum computer) that will outperform any 
classical computing device. 

Recently, we have started to realize that entanglement 
might also be related to thermodynamic properties of 
many-body systems pQ although we are far from having 
a complete understanding of this relation. 

In this short note we try to establish yet another ev- 
idence linking entanglement with thermodynamic prop- 
erties of large many body systems in the context of the 
so-called Off-Diagonal-Long- Range-Order (ODLRO), the 
concept introduced in the Ref. [2]. 

The paper is organized as follows. In the first section 
we discuss the concept of entanglement of massive par- 
ticles with vacuum. The aim of this section is to show 
that this kind of entanglement is physical, i.e., it can be 
used for quantum information processing. 

The second section discusses a particular example of 
ODLRO, namely, the one existing in the BEC phase of a 
non-interacting bosonic gas. 

Finally, in the third section we define in general the 
concept of ODLRO and argue that it is naturally linked 
to entanglement. 

I. ENTANGLEMENT OF MASSIVE PARTICLES 
WITH VACUUM 

Let us consider an infinite square well potential with 
a single particle in the first excited state inside the po- 
tential. In the second quantization formalism this can 
be written as a.| | O) , where the operator a\ creates the 
particle in the first excited state from the vacuum. If we 
now divide the well into half, creating a partition A on 
the left and a partition B on the right from the centre of 
the well, we can formally write 

a\\n)=[^=A\ + l=Bl^\n), (1) 

where the operator A\ creates a particle in the A parti- 
tion described by the sinusoidal function being the part 
of the initial wave function and the operator creates 
a particle in the B partition described by the other half 
of the initial wave function. 

The above state can be interpreted as a superposition 
of a particle being on the right (partition B) and no par- 
ticle on the left (partition A) with the particle being on 



the left (partition A) and no particle on the right (parti- 
tion B), i.e., 

lu>i«B> + iai>|iB>. 

The kets \£Ia), \^b) denote local vacuums. This kind of 
states spur lots of controversies regarding their entangle- 
ment status [3l IU [6] . We do not want to discuss this 
issue in details here but we would like to mention that 
in the case of massless particles (photons for instance) 
there seems to be no problem as the state can be used 
to entangle two distant atoms. In other words, for mass- 
less particles one can locally swap the entanglement in 
the state to two massive particles that later can be lo- 
cally manipulated and used as a resource for quantum 
information processing. The photon is destroyed in the 
process. 

It has been already shown in the Ref. [7] that one can 
extract entanglement from a single massive particle as 
well and here we briefly present slightly expanded and 
modified approach presented there. 

To this end imagine that we have two boxes each con- 
taining a single particle in the ground state of the box. 
The box A will interact only with the partition A and 
the box B with the partition B. The interaction is given 

by 

v = (a c a m\ + a\c\m a ) + (b q c b m b + bIc%m b ), 

(2) 

where the operators C A and Cb destroy the particle in 
the box A and B and the operators M A ,M B create a 
molecule in the respective box. The molecule is created 
from the merger of the particle in the box and the par- 
ticle in the well. The hamiltonian of the whole system 
consists of the free evolutions of each subsystem and the 
interaction V. We assume that the initial state is given 

by 

\i/> i )=altf A tf B \n), (3) 

i.e., there is one particle in each box and there are no 
molecules. 

In our scenario, we rapidly switch on and off the cou- 
pling (Dirac's delta coupling) V such that the whole evo- 
lution reads 

IV/) = e* v ^<) = 

{cos(g)AlC A C B + cos(g)BiC A C B + 
ism(g)C B M A + ism(g)C A M B ) (4) 
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In the Fock representation, the final state reads 

\4>f) = -^cob( S )(|1,0) + |0,1))|1, 0,1,0) + 

i»Bin( ff )|0,0)(|0 ) 1,1,0) + |1,0,1,0)), (5) 

where the ket \k, I, m, n, o,p) denotes k particles in the 
partition A, I particles in the partition B, m particles 
in the box A, n molecules in the box A, o particles in 
the box B and p molecules in the box B. Before tracing 
out the particle in the well to obtain the reduced density 
matrix of the system of particles and molecules in the 
boxes A and B let us introduce a convenient notation. 
The ket representing one particle in the box A (B) and 
no molecules in the box A (B) is denoted by | T)a(b) 
and the ket representing no particles in the box A (B) 
and one molecule in the box A (B) by | I)a(b)- I n this 
notation, the reduced density matrix after the particle in 
the well has been traced out reads 

cos( ff ) 2 | T)aa(T I® I T)sb(T 1 + 
sm(g) 2 \^) ABAB (^_\, (6) 

where |V-) = ^ (| T)a| 1)b + \ 1)a\ T)s)- The negativ- 
ity of this state is given by 

- (cos( 5 ) 2 - v/cosG?) 4 + sin( ff ) 4 ) 

and is always non-zero as long as g does not equal to the 
multiplicity of tt. For g = ^ we get maximal entangle- 
ment. The same extraction scheme works for an arbitrary 
wave function describing the particle in the box. 

Interestingly, after the interaction extracting the max- 
imal amount of entanglement the particle in the well dis- 
appears in the exact analogy with the scheme of extract- 
ing entanglement from a single photon. 

We see that in this simple scenario we have locally 
extracted entanglement carried by a single massive par- 
ticle and swapped it to a composite system of two local 
subsystems each containing massive particles of different 
kind. To verify this entanglement one must be able to 
measure the phase between the superposition of two dif- 
ferent particles as mentioned in the Ref. [7J . This can be 
done with the methods described in the Ref. [5]. 

II. BEC PHASE TRANSITION 

The presented scenario shows that the entanglement 
of a single massive particle is as physical as any other 
form of entanglement, which has been already argued 
in the Ref. [TJ. Let us know apply this results to a 
system of non-interacting spinless bosons in some three 
dimensional trapping potential. 

It is known, see the Ref. [9] for instance, that the 
essential properties of a non-interacting gas of bosons 
near the point of the BEC phase transition are contained 



in the single particle density matrix. For a canonical 
ensemble of N bosons, this density matrix reads 

k 

where the vector k labels the eigenstates of the potential 
(in the case of the three dimensional box these are the 
momenta of the particle) and (ng) is the mean number 

of bosons occupying the energy level associated with k. 
For each k the wave function \<p^) can be written as 

l^) = ^l^)l^) + ^l^)l^)' (8) 

where the ket \A^) ((A^\A^) — 1) represents the part 
of the wave function \(f)£) located in the half of the box 
(partition A) and the ket \£l B ) represents the vacuum in 
the other half of the box (partition B) , and vice versa. 

Now, we are ready to perform partial transposition on 
pi with respect to, say, partition B of the box. This par- 
tial transposition has a well defined meaning as it detects 
extractable entanglement. We get 

pI B = \Y / (\Aj:)(A j :\®\n B )(n B \ 

k 

+\il A ){Cl A \®\B i: ){B jl \) + 

^(\n A ,n B )( x \ + \x)(n A ,n B \), (9) 

where the ket \x) reads 

k 

and the norm squared of \x) is given by 

k,T 

where we have assumed that we deal with real wave func- 
tions. 

The only negative eigenvalue of p\ B is \ \J (xlx) & n d it 
can be plotted against the temperature once one knows 
the overlaps (A^\Aj), (B^\Bj). However, it is already 
clear that the plot will sharply change around the crit- 
ical temperature for condensation because of the domi- 
nance of the term near the critical temperature for 
N — > oo. This observation establishes a simple connec- 
tion between a single particle entanglement and Bose- 
Einstein condensation. In other words, one can quantify 
the BEC phase transition in terms of single-particle en- 
tanglement but this is a simple consequence of the fact 
that the single density matrix is enough to describe the 
properties of the whole system near the point of the phase 
transition. 
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We would like to stress here that the spatial entangle- 
ment and its connection to BEC have been investigated 
in the Refs. (TUl HD H2J H3] m the context of entangle- 
ment of all the particles in the gas distributed amongst 
two partitions. The conclusions in the Ref. [TJ] con- 
cerning a non-interacting gas of bosons are exactly the 
same. 

The presented results do not depend on the division 
of the potential into two parts as long as the partitions 
do not overlap. However, if the partitions A and B are 
not adjacent, i.e., there is a partition C between them, 
the reduced single-particle density matrix will contain a 
sector of vacuum in both partitions A and B, which will 
decrease the amount of entanglement. More precisely, 
the negativity in such a situation reads 



1 



(q-VmV + i 2 ), 



(12) 



where 



being the probability of finding a particle with the 
wave vector k in the partition M and 

("•£> (C) 

If the gap between the partition A and B, i.e., the size of 
the partition C, is too big q could be much larger than 
(S\6) reducing the amount of entanglement. Interestingly, 
in the absolute zero temperature the negativity reads 



N 



P?W>\A S )\B S ) 



(C) 



\ 



(-4) (B) 

■ go Po 

(C) (C) 

Po Po 



where p M ^ is the probability of finding a particle in the 
ground state in the partition M. This expression is never 
zero as long as there is a chance of finding a particle in 
the partitions A and B. 



III. RELATION OF ODLRO TO 
ENTANGLEMENT 

Let us know generalize the above results to demon- 
strate a link between the concept of ODLRO and en- 
tanglement. Imagine a system of N interacting parti- 
cles (bosons or fermions). According to the Ref. [2] one 
can identify the different thermodynamical phases of the 
system by looking at the spectral properties of the re- 
duced ro-body density matrices p n . For what follows next 
we set the normalization of the density operators p n to 

The first reduced density matrix to investigate is p\. 
The ODLRO in this matrix is defined as the existence 



of an eigenvalue of the order of N. This automatically 
means that the other eigenvalues are small and p\ can be 
written as 



Pl = aN\^\+p' x 



(13) 



where p\ is a small positive operator and a is some pos- 
itive number. As we did before, we can quantify the 
single-particle entanglement by the negativity of the par- 
tially transposed matrix p\ with respect to some partition 
of the trapping potential into two disjoined parts. We can 
expect, as in the case of the BEC discussed earlier, that 
the entanglement will sharply change around the point 
where the spectrum p\ becomes concentrated around the 
dominating eigenvalue. 

The intuitive understanding of this behaviour of en- 
tanglement is the following one. As was shown in the 
Ref. [2] the existence of the eigenvalue of the order of TV 
implies that 



aJV 
~V~ 



(14) 



for \r — s\ — > oo, where V is the volume occupied by the 
system of N bosons. The above equation states that the 
amplitude of the probability of the process describing 
hopping of a particle from the point r to the point s 
becomes substantial. This, in turn, implies that one can 
expect to have a single-particle entanglement between 
these two points in space; one has the superposition of 
the particle at point r and no particle at point s with the 
particle at point s and no particle at point r. Therefore, 
whenever there is ODLRO in the single-particle reduced 
density matrix there is a substantial amount of the single- 
particle entanglement. 

One can also talk about ODLRO present in the two- 
body reduced density matrix without ODLRO in the one- 
body density matrix j2j. Its presence is detected by the 
existence of an eigenvalue of the order of N with all the 
other ones being of the lower order (a good physical ex- 
ample of this kind is low temperature superconductivity, 
the BCS state). In this case we can write 



p 2 = aN\cj ) ){4 > \ + p' 2l 



(15) 



where \<j>) is some two particle wave function and, as be- 
fore, the p' 2 is a small positive operator. In the position 
representation this condition means [2J that only the di- 
agonal elements of ^(^l^l^i^a) and off diagonal ones of 
the form p2(rr\ss) are non-zero, even for |r — s\ — > oo. 
Strictly speaking, this happens for r and s vectors micro- 
scopically close to each other; they do not have to point 
to the same point in space. Physically, this condition 
means that the amplitude of probability of two particles 
hopping from one location in space to another one dom- 
inates all the other processes. As the consequence, we 
expect to have entanglement generated between two ar- 
bitrary partitions. This entanglement will be of the form 



)a\0)b + \0)a\(I> 1 



(16) 
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where the ket \4>a) represents the part of the two-particle 
ket \<j)) describing two particles in the partition A (the 
same applies to the partition B). Additionally, this en- 
tanglement becomes substantial only in the regime where 
ODLRO in p2 exists, for instance in the superconducting 
states. 

In this short note we have presented a straightforward 
connection between the concept of ODLRO for a system 
of N interacting bosons or fermions and entanglement. It 
seems that ODLRO implies entanglement of massive par- 
ticles with vacuum (at least in the cases discussed in this 
note) but we do not know if the presence of this kind 
of entanglement generates ODLRO. It would be inter- 



esting to investigate higher order ODLRO, i.e., the one 
present in n-body reduced density matrices and its rela- 
tion to entanglement. This will be the subject of further 
research. 
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